RPA and Rad27 limit templated and inverted insertions at DNA breaks.

Formation of templated insertions at DNA double-strand breaks (DSBs) is very common in cancer cells. The mechanisms and enzymes regulating these events are largely unknown. Here, we investigated templated insertions in yeast at DSBs using amplicon sequencing across a repaired locus. We document very short (most ∼5-34 bp), templated inverted duplications at DSBs. They are generated through a foldback mechanism that utilizes microhomologies adjacent to the DSB. Enzymatic requirements suggest a hybrid mechanism wherein one end requires Polδ-mediated synthesis while the other end is captured by nonhomologous end joining (NHEJ). This process is exacerbated in mutants with low levels or mutated RPA ( rtt105 Δ; rfa1 -t33) or extensive resection mutant ( sgs1 Δ exo1 Δ). Templated insertions from various distant genomic locations also increase in these mutants as well as in rad27 Δ and originate from fragile regions of the genome. Among complex insertions, common events are insertions of two sequences, originating from the same locus and with inverted orientation. We propose that these inversions are also formed by microhomology-mediated template switching. Taken together, we propose that a shortage of RPA typical in cancer cells is one possible factor stimulating the formation of templated insertions.

Yeast EndoG prevents genome instability by degrading cytoplasmic DNA.

In metazoans release of mitochondrial DNA or retrotransposon cDNA to cytoplasm can cause sterile inflammation and disease 1. Cytoplasmic nucleases degrade these DNA species to limit inflammation 2,3. It remains unknown whether degradation these DNA also prevents nuclear genome instability. To address this question, we decided to identify the nuclease regulating transfer of these cytoplasmic DNA species to the nucleus. We used an amplicon sequencing-based method in yeast enabling analysis of millions of DSB repair products. Nuclear mtDNA (NUMTs) and retrotransposon cDNA insertions increase dramatically in nondividing stationary phase cells. Yeast EndoG (Nuc1) nuclease limits insertions of cDNA and transfer of very long mtDNA (>10 kb) that forms unstable circles or rarely insert in the genome, but it promotes formation of short NUMTs (~45-200 bp). Nuc1 also regulates transfer of cytoplasmic DNA to nucleus in aging or during meiosis. We propose that Nuc1 preserves genome stability by degrading retrotransposon cDNA and long mtDNA, while short NUMTs can originate from incompletely degraded mtDNA. This work suggests that nucleases eliminating cytoplasmic DNA play a role in preserving genome stability.

Drugging evolution of antibiotic resistance at a regulatory network hub.

Evolution of antibiotic resistance is a world health crisis, fueled by new mutations. Drugs to slow mutagenesis could, as cotherapies, prolong the shelf-life of antibiotics, yet evolution-slowing drugs and drug targets have been underexplored and ineffective. Here, we used a network-based strategy to identify drugs that block hubs of fluoroquinolone antibiotic-induced mutagenesis. We identify a U.S. Food and Drug Administration- and European Medicines Agency-approved drug, dequalinium chloride (DEQ), that inhibits activation of the Escherichia coli general stress response, which promotes ciprofloxacin-induced (stress-induced) mutagenic DNA break repair. We uncover the step in the pathway inhibited: activation of the upstream "stringent" starvation stress response, and find that DEQ slows evolution without favoring proliferation of DEQ-resistant mutants. Furthermore, we demonstrate stress-induced mutagenesis during mouse infections and its inhibition by DEQ. Our work provides a proof-of-concept strategy for drugs to slow evolution in bacteria and generally.

ppGpp and RNA-polymerase backtracking guide antibiotic-induced mutable gambler cells.

Antibiotic resistance is a global health threat and often results from new mutations. Antibiotics can induce mutations via mechanisms activated by stress responses, which both reveal environmental cues of mutagenesis and are weak links in mutagenesis networks. Network inhibition could slow the evolution of resistance during antibiotic therapies. Despite its pivotal importance, few identities and fewer functions of stress responses in mutagenesis are clear. Here, we identify the Escherichia coli stringent starvation response in fluoroquinolone-antibiotic ciprofloxacin-induced mutagenesis. Binding of response-activator ppGpp to RNA polymerase (RNAP) at two sites leads to an antibiotic-induced mutable gambler-cell subpopulation. Each activates a stress response required for mutagenic DNA-break repair: surprisingly, ppGpp-site-1-RNAP triggers the DNA-damage response, and ppGpp-site-2-RNAP induces σS-response activity. We propose that RNAP regulates DNA-damage processing in transcribed regions. The data demonstrate a critical node in ciprofloxacin-induced mutagenesis, imply RNAP-regulation of DNA-break repair, and identify promising targets for resistance-resisting drugs.

Yeast EndoG prevents genome instability by degrading cytoplasmic DNA.

In metazoans release of mitochondrial DNA or retrotransposon cDNA to cytoplasm can cause sterile inflammation and disease. Cytoplasmic nucleases degrade these DNA species to limit inflammation. It remains unknown whether degradation these DNA also prevents nuclear genome instability. To address this question, we decided to identify the nuclease regulating transfer of these cytoplasmic DNA species to the nucleus. We used an amplicon sequencing-based method in yeast enabling analysis of millions of DSB repair products. Nu clear mt DNA (NUMTs) and retrotransposon cDNA insertions increase dramatically in nondividing stationary phase cells. Yeast EndoG (Nuc1) nuclease limits insertions of cDNA and transfer of very long mtDNA (>10 kb) that forms unstable circles or rarely insert in the genome, but it promotes formation of short NUMTs (∼45-200 bp). Nuc1 also regulates transfer of cytoplasmic DNA to nucleus in aging or during meiosis. We propose that Nuc1 preserves genome stability by degrading retrotransposon cDNA and long mtDNA, while short NUMTs can originate from incompletely degraded mtDNA. This work suggests that nucleases eliminating cytoplasmic DNA play a role in preserving genome stability.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation.

BACKGROUND: The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases. METHODS: Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype. RESULTS: In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations. CONCLUSIONS: Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Stress-Induced Mutagenesis, Gambler Cells, and Stealth Targeting Antibiotic-Induced Evolution.

Mechanisms of evolution and evolution of antibiotic resistance are both fundamental and world health problems. Stress-induced mutagenesis defines mechanisms of mutagenesis upregulated by stress responses, which drive adaptation when cells are maladapted to their environments-when stressed. Work in mutagenesis induced by antibiotics had produced tantalizing clues but not coherent mechanisms. We review recent advances in antibiotic-induced mutagenesis that integrate how reactive oxygen species (ROS), the SOS and general stress responses, and multichromosome cells orchestrate a stress response-induced switch from high-fidelity to mutagenic repair of DNA breaks. Moreover, while sibling cells stay stable, a mutable "gambler" cell subpopulation is induced by differentially generated ROS, which signal the general stress response. We discuss other evolvable subpopulations and consider diverse evolution-promoting molecules as potential targets for drugs to slow evolution of antibiotic resistance, cross-resistance, and immune evasion. An FDA-approved drug exemplifies "stealth" evolution-slowing drugs that avoid selecting resistance to themselves or antibiotics.

Two mechanisms of chromosome fragility at replication-termination sites in bacteria.

Chromosomal fragile sites are implicated in promoting genome instability, which drives cancers and neurological diseases. Yet, the causes and mechanisms of chromosome fragility remain speculative. Here, we identify three spontaneous fragile sites in the Escherichia coli genome and define their DNA damage and repair intermediates at high resolution. We find that all three sites, all in the region of replication termination, display recurrent four-way DNA or Holliday junctions (HJs) and recurrent DNA breaks. Homology-directed double-strand break repair generates the recurrent HJs at all of these sites; however, distinct mechanisms of DNA breakage are implicated: replication fork collapse at natural replication barriers and, unexpectedly, frequent shearing of unsegregated sister chromosomes at cell division. We propose that mechanisms such as both of these may occur ubiquitously, including in humans, and may constitute some of the earliest events that underlie somatic cell mosaicism, cancers, and other diseases of genome instability.

Gamblers: An Antibiotic-Induced Evolvable Cell Subpopulation Differentiated by Reactive-Oxygen-Induced General Stress Response.

Antibiotics can induce mutations that cause antibiotic resistance. Yet, despite their importance, mechanisms of antibiotic-promoted mutagenesis remain elusive. We report that the fluoroquinolone antibiotic ciprofloxacin (cipro) induces mutations by triggering transient differentiation of a mutant-generating cell subpopulation, using reactive oxygen species (ROS). Cipro-induced DNA breaks activate the Escherichia coli SOS DNA-damage response and error-prone DNA polymerases in all cells. However, mutagenesis is limited to a cell subpopulation in which electron transfer together with SOS induce ROS, which activate the sigma-S (σS) general-stress response, which allows mutagenic DNA-break repair. When sorted, this small σS-response-"on" subpopulation produces most antibiotic cross-resistant mutants. A U.S. Food and Drug Administration (FDA)-approved drug prevents σS induction, specifically inhibiting antibiotic-promoted mutagenesis. Further, SOS-inhibited cell division, which causes multi-chromosome cells, promotes mutagenesis. The data support a model in which within-cell chromosome cooperation together with development of a "gambler" cell subpopulation promote resistance evolution without risking most cells.

Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2.

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

Bacteria-to-Human Protein Networks Reveal Origins of Endogenous DNA Damage.

DNA damage provokes mutations and cancer and results from external carcinogens or endogenous cellular processes. However, the intrinsic instigators of endogenous DNA damage are poorly understood. Here, we identify proteins that promote endogenous DNA damage when overproduced: the DNA "damage-up" proteins (DDPs). We discover a large network of DDPs in Escherichia coli and deconvolute them into six function clusters, demonstrating DDP mechanisms in three: reactive oxygen increase by transmembrane transporters, chromosome loss by replisome binding, and replication stalling by transcription factors. Their 284 human homologs are over-represented among known cancer drivers, and their RNAs in tumors predict heavy mutagenesis and a poor prognosis. Half of the tested human homologs promote DNA damage and mutation when overproduced in human cells, with DNA damage-elevating mechanisms like those in E. coli. Our work identifies networks of DDPs that provoke endogenous DNA damage and may reveal DNA damage-associated functions of many human known and newly implicated cancer-promoting proteins.

Fluorescent fusions of the N protein of phage Mu label DNA damage in living cells.

The N protein of phage Mu was indicated from studies in Escherichia coli to hold linear Mu chromosomes in a circular conformation by non-covalent association, and thus suggested potentially to bind DNA double-stranded ends. Because of its role in association with linear Mu DNA, we tested whether fluorescent-protein fusions to N might provide a useful tool for labeling DNA damage including double-strand break (DSB) ends in single cells. We compared N-GFP with a biochemically well documented DSB-end binding protein, the Gam protein of phage Mu, also fused to GFP. We find that N-GFP produced in live E. coli forms foci in response to DNA damage induced by radiomimetic drug phleomycin, indicating that it labels damaged DNA. N-GFP also labels specific DSBs created enzymatically by I-SceI double-strand endonuclease, and by X-rays, with the numbers of foci corresponding with the numbers of DSBs generated, indicating DSB labeling. However, whereas N-GFP forms about half as many foci as GamGFP with phleomycin, its labeling of I-SceI- and X-ray-induced DSBs is far less efficient than that of GamGFP. The data imply that N-GFP binds and labels DNA damage including DSBs, but may additionally label phleomycin-induced non-DSB damage, with which DSB-specific GamGFP does not interact. The data indicate that N-GFP labels DNA damage, and may be useful for general, not DSB-specific, DNA-damage detection.

Oxygen and RNA in stress-induced mutation.

Mechanisms of mutation upregulated by stress responses have been described in several organisms from bacteria to human. These mechanisms might accelerate genetic change specifically when cells are maladapted to their environment. Stress-induced mutation mechanisms differ in their genetic requirements from mutation in growing cells, occurring by different mechanisms in different assay systems, but having in common a requirement for the induction of stress-responses. Here, we review progress in two areas relevant to stress-response-dependent mutagenic DNA break repair mechanisms in Escherichia coli. First, we review evidence that relates mutation to transcription. This connection might allow mutagenesis in transcribed regions, including those relevant to any stress being experienced, opening the possibility that mutations could be targeted to regions where mutation might be advantageous under conditions of a specific stress. We review the mechanisms by which replication initiated by transcription can lead to mutation. Second, we review recent findings that, although stress-induced mutation does not require exogenous DNA-damaging agents, it does require the presence of damaged bases in DNA. For starved E. coli, endogenous oxygen radicals cause these altered bases. We postulate that damaged bases stall the replisome, which, we suggest, is required for DNA-polymerase exchange, allowing the action of low-fidelity DNA polymerases that promote mutation.

Persistent damaged bases in DNA allow mutagenic break repair in Escherichia coli.

Bacteria, yeast and human cancer cells possess mechanisms of mutagenesis upregulated by stress responses. Stress-inducible mutagenesis potentially accelerates adaptation, and may provide important models for mutagenesis that drives cancers, host pathogen interactions, antibiotic resistance and possibly much of evolution generally. In Escherichia coli repair of double-strand breaks (DSBs) becomes mutagenic, using low-fidelity DNA polymerases under the control of the SOS DNA-damage response and RpoS general stress response, which upregulate and allow the action of error-prone DNA polymerases IV (DinB), II and V to make mutations during repair. Pol IV is implied to compete with and replace high-fidelity DNA polymerases at the DSB-repair replisome, causing mutagenesis. We report that up-regulated Pol IV is not sufficient for mutagenic break repair (MBR); damaged bases in the DNA are also required, and that in starvation-stressed cells, these are caused by reactive-oxygen species (ROS). First, MBR is reduced by either ROS-scavenging agents or constitutive activation of oxidative-damage responses, both of which reduce cellular ROS levels. The ROS promote MBR other than by causing DSBs, saturating mismatch repair, oxidizing proteins, or inducing the SOS response or the general stress response. We find that ROS drive MBR through oxidized guanines (8-oxo-dG) in DNA, in that overproduction of a glycosylase that removes 8-oxo-dG from DNA prevents MBR. Further, other damaged DNA bases can substitute for 8-oxo-dG because ROS-scavenged cells resume MBR if either DNA pyrimidine dimers or alkylated bases are induced. We hypothesize that damaged bases in DNA pause the replisome and allow the critical switch from high fidelity to error-prone DNA polymerases in the DSB-repair replisome, thus allowing MBR. The data imply that in addition to the indirect stress-response controlled switch to MBR, a direct cis-acting switch to MBR occurs independently of DNA breakage, caused by ROS oxidation of DNA potentially regulated by ROS regulators.

Stress-Induced Mutagenesis: Implications in Cancer and Drug Resistance.

Genomic instability underlies many cancers and generates genetic variation that drives cancer initiation, progression, and therapy resistance. In contrast with classical assumptions that mutations occur purely stochastically at constant, gradual rates, microbes, plants, flies, and human cancer cells possess mechanisms of mutagenesis that are upregulated by stress responses. These generate transient, genetic-diversity bursts that can propel evolution, specifically when cells are poorly adapted to their environments-that is, when stressed. We review molecular mechanisms of stress-response-dependent (stress-induced) mutagenesis that occur from bacteria to cancer, and are activated by starvation, drugs, hypoxia, and other stressors. We discuss mutagenic DNA break repair in Escherichia coli as a model for mechanisms in cancers. The temporal regulation of mutagenesis by stress responses and spatial restriction in genomes are common themes across the tree of life. Both can accelerate evolution, including the evolution of cancers. We discuss possible anti-evolvability drugs, aimed at targeting mutagenesis and other variation generators, that could be used to delay the evolution of cancer progression and therapy resistance.

An ultra-dense library resource for rapid deconvolution of mutations that cause phenotypes in Escherichia coli.

With the wide availability of whole-genome sequencing (WGS), genetic mapping has become the rate-limiting step, inhibiting unbiased forward genetics in even the most tractable model organisms. We introduce a rapid deconvolution resource and method for untagged causative mutations after mutagenesis, screens, and WGS in Escherichia coli. We created Deconvoluter-ordered libraries with selectable insertions every 50 kb in the E. coli genome. The Deconvoluter method uses these for replacement of untagged mutations in the genome using a phage-P1-based gene-replacement strategy. We validate the Deconvoluter resource by deconvolution of 17 of 17 phenotype-altering mutations from a screen of N-ethyl-N-nitrosourea-induced mutants. The Deconvoluter resource permits rapid unbiased screens and gene/function identification and will enable exploration of functions of essential genes and undiscovered genes/sites/alleles not represented in existing deletion collections. This resource for unbiased forward-genetic screens with mapping-by-sequencing ('forward genomics') demonstrates a strategy that could similarly enable rapid screens in many other microbes.

Holliday junction trap shows how cells use recombination and a junction-guardian role of RecQ helicase.

DNA repair by homologous recombination (HR) underpins cell survival and fuels genome instability, cancer, and evolution. However, the main kinds and sources of DNA damage repaired by HR in somatic cells and the roles of important HR proteins remain elusive. We present engineered proteins that trap, map, and quantify Holliday junctions (HJs), a central DNA intermediate in HR, based on catalytically deficient mutant RuvC protein of Escherichia coli. We use RuvCDefGFP (RDG) to map genomic footprints of HR at defined DNA breaks in E. coli and demonstrate genome-scale directionality of double-strand break (DSB) repair along the chromosome. Unexpectedly, most spontaneous HR-HJ foci are instigated, not by DSBs, but rather by single-stranded DNA damage generated by replication. We show that RecQ, the E. coli ortholog of five human cancer proteins, nonredundantly promotes HR-HJ formation in single cells and, in a novel junction-guardian role, also prevents apparent non-HR-HJs promoted by RecA overproduction. We propose that one or more human RecQ orthologs may act similarly in human cancers overexpressing the RecA ortholog RAD51 and find that cancer genome expression data implicate the orthologs BLM and RECQL4 in conjunction with EME1 and GEN1 as probable HJ reducers in such cancers. Our results support RecA-overproducing E. coli as a model of the many human tumors with up-regulated RAD51 and provide the first glimpses of important, previously elusive reaction intermediates in DNA replication and repair in single living cells.

Roles of Nucleoid-Associated Proteins in Stress-Induced Mutagenic Break Repair in Starving Escherichia coli.

The mutagenicity of DNA double-strand break repair in Escherichia coli is controlled by DNA-damage (SOS) and general (RpoS) stress responses, which let error-prone DNA polymerases participate, potentially accelerating evolution during stress. Either base substitutions and indels or genome rearrangements result. Here we discovered that most small basic proteins that compact the genome, nucleoid-associated proteins (NAPs), promote or inhibit mutagenic break repair (MBR) via different routes. Of 15 NAPs, H-NS, Fis, CspE, and CbpA were required for MBR; Dps inhibited MBR; StpA and Hha did neither; and five others were characterized previously. Three essential genes were not tested. Using multiple tests, we found the following: First, Dps, which reduces reactive oxygen species (ROS), inhibited MBR, implicating ROS in MBR. Second, CbpA promoted F' plasmid maintenance, allowing MBR to be measured in an F'-based assay. Third, Fis was required for activation of the SOS DNA-damage response and could be substituted in MBR by SOS-induced levels of DinB error-prone DNA polymerase. Thus, Fis promoted MBR by allowing SOS activation. Fourth, H-NS represses ROS detoxifier sodB and was substituted in MBR by deletion of sodB, which was not otherwise mutagenic. We conclude that normal ROS levels promote MBR and that H-NS promotes MBR by maintaining ROS. CspE positively regulates RpoS, which is required for MBR. Four of five previously characterized NAPs promoted stress responses that enhance MBR. Hence, most NAPs affect MBR, the majority via regulatory functions. The data show that a total of six NAPs promote MBR by regulating stress responses, indicating the importance of nucleoid structure and function to the regulation of MBR and of coupling mutagenesis to stress, creating genetic diversity responsively.

Absence of heterozygosity due to template switching during replicative rearrangements.

We investigated complex genomic rearrangements (CGRs) consisting of triplication copy-number variants (CNVs) that were accompanied by extended regions of copy-number-neutral absence of heterozygosity (AOH) in subjects with multiple congenital abnormalities. Molecular analyses provided observational evidence that in humans, post-zygotically generated CGRs can lead to regional uniparental disomy (UPD) due to template switches between homologs versus sister chromatids by using microhomology to prime DNA replication-a prediction of the replicative repair model, MMBIR. Our findings suggest that replication-based mechanisms might underlie the formation of diverse types of genomic alterations (CGRs and AOH) implicated in constitutional disorders.

Engineered proteins detect spontaneous DNA breakage in human and bacterial cells.

Spontaneous DNA breaks instigate genomic changes that fuel cancer and evolution, yet direct quantification of double-strand breaks (DSBs) has been limited. Predominant sources of spontaneous DSBs remain elusive. We report synthetic technology for quantifying DSBs using fluorescent-protein fusions of double-strand DNA end-binding protein, Gam of bacteriophage Mu. In Escherichia coli GamGFP forms foci at chromosomal DSBs and pinpoints their subgenomic locations. Spontaneous DSBs occur mostly one per cell, and correspond with generations, supporting replicative models for spontaneous breakage, and providing the first true breakage rates. In mammalian cells GamGFP-labels laser-induced DSBs antagonized by end-binding protein Ku; co-localizes incompletely with DSB marker 53BP1 suggesting superior DSB-specificity; blocks resection; and demonstrates DNA breakage via APOBEC3A cytosine deaminase. We demonstrate directly that some spontaneous DSBs occur outside of S phase. The data illuminate spontaneous DNA breakage in E. coli and human cells and illustrate the versatility of fluorescent-Gam for interrogation of DSBs in living cells. DOI:http://dx.doi.org/10.7554/eLife.01222.001.